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ABSTRACT 

Context. Contrary to theoretical expectation, surprisingly low concentrations of molecular oxygen, 2 , have been found in the inter- 
stellar medium by means of orbiting telescopes. 

Aims. Observations of the (Nj = \\ - 1 ) ground state transition of 2 with the Odin satellite resulted in a > 5<r detection toward the 
dense core pOphA. At the frequency of the line, 119 GHz, the Odin telescope has a beam width of 10', larger than the size of the 
dense core, so that the precise nature of the emitting source and its exact location and extent are unknown. The current investigation 
is intended to remedy this. 

Methods. Telluric absorption makes ground based 2 observations essentially impossible and observations had to be done from space. 
Millimetre-wave telescopes on space platforms were necessarily small, which resulted in large, several arcminutes wide, beam pat- 
terns. Although the Earth's atmosphere is entirely opaque to low-lying 2 transitions, it allows ground based observations of the much 
rarer 16 18 in favourable conditions and at much higher angular resolution with larger telescopes. In addition, p Oph A exhibits both 
multiple radial velocity systems and considerable velocity gradients. Extensive mapping of the region in the proxy C 18 (/ = 3 - 2) 
line can be expected to help identify the 2 source on the basis of its line shape and Doppler velocity. Line opacities were determined 
from observations of optically thin 13 C 18 (/ = 3 - 2) at selected positions. 

Results. During several observing periods, two C 18 intensity maxima in pOph A were searched for in the 16 18 (2i - 0i) line at 
234 GHz with the 12 m APEX telescope. These positions are associated with peaks in the mm-continuum emission from dust. Our 
observations resulted in an upper limit on the integrated 18 intensity of fT^dv < 0.01 Kkm s _1 (3cr) into the 26'/5 beam. 
Conclusions. Examining the evidence, which is based primarily on observations in lines of 18 and C 18 0, leads us to conclude that 
the source of observed 2 emission is most likely confined to the central regions of the p Oph A core. In this limited area, implied 2 
abundances could thus be higher than previously reported, by up to two orders of magnitude. 
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1. Introduction 

Oxy gen is the most abundant of the astronomical metals (e.g., 
|Asplund et al. 2009 1 and references therein). Consequently, in 
its molecular form, it was also expected to be very abundant in 
the UV-shielded regions inside molecular clouds (e.g., Black & 
Simth||1984[|Bergin et aL||2000| piarnley et aL||2001[ |Roberts 



& Herbst|266 2; Spaan s & van D ishoeck p001[|Viti et al.|2001[ 
Willacy et al.|2002[|Quan et al.|2008| ) and to contribute signifi 



cantly to the cooling, hence the energy balance, of dense clouds 
( [Goldsmith & Langer|1978 ). 

Because of the high O2 content in the Earth's atmosphere, 
astronomical O2 sources cannot be observed from the ground. 
Dedicated space mission^] came into operation near the be- 
ginning of the new millenium. Their unsuccessful searches 
( [Goldsmith eTaL][2000l |Pagani et aL]|2003] ) were highly disap- 
pointing and it was hard to understand that, in the interstellar 
medium (ISM), O2 is an elusive species (see references cited 
above). 



Send offprint requests to: R. Liseau 

* Based on observations with APEX, Llano Chajnantor, Chile. 

1 SWAS in 1998, see http : //cf a-www . harvard . edu/swas/ ., and 
Odin in 2001, see http : //www . snsb . se/eng_odin_intro . shtml 



Eventually, after more than 20 days of Odin-observing dur- 
ing three different runs, came a real break-throug h: for the very 
first time, O2 was finally detected in the ISM ( [Larsson et al. 
200 7). The O2 emitting object, p Oph A, is a dense clump ( [Loren 
et al.|l990 ) in a region of active star formation (L 1688). On the 



basis of theoretical model calculations, the detection on the basis 
of theoretical model calculationsof this kind of source had ear- 



lier b een predicted by Black & Smith] ( [1984] ) and |Marechal et al. 
(1997a), where the latter authors made their specific prediction 
with regard to Odin. 

Odin carries a 1 . 1 m telescope which is designed for obser- 
vations in the submillimetre regime, between roughly 480 and 
580 GHz (0.5 - 0.6 mm). However, the O2 discovery was made 
with a dedicated 119 GHz (2.5 mm) receiver aboard Odin, fix- 
tuned to the frequency of the ground state O2 (Nj = l\ - lo) 
transition at 118 750.343 MHz. At this frequency, the telescope 
beam size is 10', larger than the angular dimension of the dense 
p Oph A core, which is about 4' (FWHM of devonvolved CS 
core, |Liseau et al.|1995| ). 

It follows that the true O2 source is likely under-resolved, 
the consequence of which directly affects estimates of the abun- 
dance of O2, i.e. Af(02)/Af(H 2 ): depending on the adopted model, 
the Odin observations imply an abun dance which is cu rrently 
uncertain by two orders of magnitude ( Liseau et al.|2005| ). 
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In Fig. 2 of Larss on et al. ( 2007 ), the Odin-02 line is com- 
pared to transitions of other molecular species in pOphA. 
Whereas lines of H 2 and CO are optically very thick over large 
parts of the cloud and have self-absorbed profiles, the optically 
thin O2 line displays a simple, Gaussian shape. This line shape 
is similar to that of a carbon recombination line, displayed at the 
top of the figure and which most likely originates in the p Oph- 
PDR. If also the main source of O2 emission, the abundance 
would indeed be very low. 

However, the O2 line shape is also similar to that of the C 18 O 
(3-2) line, also shown in the figure. This suggests that the C 18 
line can be used as a tracer of the molecular oxygen emission 
and we set out to map the 10' Odin beam in the (3-2) transi- 
tion of C 18 with the APEX beam of size 19". It was expected 
that a detailed comparison of the line centre velocity with that of 
the O2 line would help to narrow down the exact location of the 
O2 emission, since two distinct velocity components are known 
to be present in p Oph A. This information is needed to under- 
stand, where, i.e. in what physical conditions, the majority of 
the O2 molecules is excited: in the cold and dense dark cores ( |Di| 
Fran cesco et al.|2004]), in the extende d warm Photon Dominated 
Region (PDR; Hollenbach et al.|2009|) or in the hot sho cked gas 
of the outflow from VLA 1623 ( |Liseau & Justtanont|2 009)? With 
the C 18 proxy for O2 emission, probable emission regions were 
identified, which were then observed for 16 18 (2\ - 0i). 

There exists earlier work for this line and the p Oph cloud. 
Goldsmith et al.| ( |1985| ) observed p Oph A in the same transi- 



tion and with comparable beam size (26"), albeit at an offset 
11"E and 61" N relative to the position of SM IN. They ob- 
tained T*< 120mK (l<x) over 0.34 km s" 1 . At similar chan- 
nel resol ution (0.32 km s" 1 ) and towar d essentially the same 
position, Liszt & Vanden Bout ( 1985| ) obtained an rms-noise 
value of T* < 17.5 mK with the 12 m NRAO telescope (34"). 
These papers also present energy level diagrams. Observations 
made with the 10 m telescope of the Caltech Submillimeter 
Observatory (CSO) in July 1991 and for the position 16 h 23 m 25 s , 
-24°15 / 49" (B 1950) H resulted in an rms noise temperature of 
16 mK in a 0.25 km s~^ velocity bin and of 12 mK after binning 
to 0.50 km s" 1 (E. van Dishoeck, J. Keene & T. Phillips, private 
communication). 

The derivation of molecular abundances requires knowledge 
of the H2 column density. One of the widely exploited tech- 
niques to estimate A/^Fb) is to use observations of C 18 0, the 
transitions of which in many cases can be shown to be optically 
thin. We discovered, however, that in the dense core regions of 
p Oph A, this not to be the case everywhere and that appropriate 
opacity corrections using the 13 C 18 O line needed to be made. 

This paper is organised as follows: In Sect. 2, our APEX ob- 
servations of the p Oph A cloud in transitions of 18 0, C 18 and 
13 C 18 are described. Sect. 3 presents our results, which are dis- 
cussed in Sect. 4. Finally, in Sect. 5 our main conclusions are 
briefly summarised. 



2. Observations and data reductions 

All observations have been made with the SIS receivers and 
spectrometers at the Atacama Pathfinder Experiment (APEX). 
The 12 m APEX telescope is located at an altitude of about 
5100 m on the Llano de Chajnantor in northern Chil^] The tele- 
scope pointing is accurate to 3" (rms). 



The Fast Fourier Transform Spectrometer (FFTS) was con- 
figured to have 8192 channels, which over a bandwidth of 1 GHz 
provides a resolution of 122 kHz, corresponding to 0.16kms _1 
and 0.1 1 km s" 1 at 234 GHz and 329 GHz, respectively. As fron- 
tends for these frequencies, we used APEX 1 of the Swedish 
Heterodyn e Facility Instrument (SHFI, |Vassilev et al. |2 008 ) and 
APEX 2A flRisacher et al.|2006) . 



2.1. 18 observations 

The data have been collected during three different observing 
runs in 2008 and 2009. The frequency of the (2\ - 0Q line c an 
be derived from the data given by Steinbach & Gordy|( [T975] ) as 
233946.179 MHz. At 234 GHz, the APEX beam has a half power 
beam width HPBW=26 / /5 and the main beam efficiency is 77 m b = 
0.75. The telescope was pointed toward RA=16 h 26 m 27 s 2 and 
Dec=-24°23'34" (J2000), a position which was initially chosen 
on the basis of, as it turned out, insufficiently sampled data (see 
Sect. 3.2). In addition, the strongest peak of doubly deuterated 
formaldhyde emission in the p Oph A core (P. Bergman et al., in 
preparation J] which is situated 30" south of these coordinates, 
was also observed. These positions are close to the location of 
intense mm-dust-emission (cf. Fig.|3]), i.e. the dense core SM 1 
( Motte et al. 1998| ). For the primary position, the total on-source 
integration time was 4.9 hours and the average system tempera- 



ture was r sys ~ 220 K, whereas for the -30" -position, these val- 
ues were 6.5 hours and 210 K, respectively. 



2.2. C 18 and 13 C 18 observations 

The observations were collected during two observing runs in 
2006 and 2007 at the APEX telescope. The observing mode 
was position switched raster mapping and the data were sam- 
pled according to the Nyquist criterion on a rectangular 10" grid, 
aligned with the equatorial coordinate system (200" x 200"). 
At 329 GHz, the HPBW=19" and the average system temper- 
ature was r sys =200K. The efficiencies were 7/ mb = 0.73 and 
?7Moon = 0-85 for point source and extended source calibrations, 
respectively. 

In addition, an extended raster map of the outer regions of 
pOphA was obtained on a coarser grid with 20" (full beam) 
spacings. The entire region observed is thus as large as AaxAS = 
10' x5'. 

The origin of the map is the same as that of the Odin ob- 
servations, i.e. the (0, 0) position is at RA=16 h 26 m 24?6 and 
Dec=-24°23'54" (J20 00). The same refere nce position as for 
the Odin observations ( [Larsson et al.||2007| ), viz. 15' N relative 
to the map centre, was used here for calibration purposes. In 
addition to the C 18 map, five positions were also observed in 
the (3-2) transition of the even rarer isotope 13 C 18 (Table[2]). 
Kla pper et al.| ( |2003| ) provide lab-frequencies for the (3-2) ro- 
tational transition of C 18 and 13 C 18 0, i.e., 329 330.552 MHz 
and 314119.660 MHz, respectively, and where the latter is a 
weighted mean value, with the 13 C hyperfine structure being ig- 
nored. 



2 This corresponds to 16 h 26 m 26 s 4, -24°22'33" in J2000 coordinates 
and is at (+25", +80") relative to the origin of the C ls O map (Fig.[3|. 

3 |http : //www . apex-tel escope . org/ 



4 The 234 GHz spectra admitted also lines of deuterated formalde- 
hyde. Mapping observations revealed this peak position. 
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233.85 233.9 233.95 234 

Frequency [GHz] 

Fig. 1. The central part of the 1 GHz wide APEX spectrum 
centered on the frequency of the 16 18 0(2i - Oi) transi- 
tion, 233.946179 GHz, and obtained toward RA=16 h 26 m 27 s 2 
and Dec=-24°24'04" (J2000) in pOphA. The sampling is in 
122 kHz wide channels (6v = 0.16 km s" 1 ). 

Table 1. C 18 0-peaks of integrated intensity, fT^dv 



C 18 0-peak 


Offset (arcsec) 


R.A. (J2000) 


Dec. (J2000) 


Other ID 


PI 


+35, +50 


16:26:27.2 


-24:23:04 


Nl 


P2 


+45, +28 


16:26:27.9 


-24:23:26 


N5, SM1N 


P3 


+45, -13 


16:26:27.9 


-24:24:07 


N4, SMI 


P4 


+48, +75 


16:26:28.1 


-24:22:39 


16264-2422b 



Note to the T able: N-source s from |Di Francesco et al.| (2004) , SM- 
objects from |Motte et al.| (1998) , RA-Dec labelled source from 
IJohnstone et al.| ( |2000) . 



3. Results 

3.1. 18 

The 16 18 (2i - 0i) line was not detected toward any of the ob- 
served positions. Toward the position associated with P2 (see 
Fig.[3] and Table[j}, the noise level is r rms =6.5mK (lcr) in a 
0.62 km s" 1 bin. The result is similar for the observation of the 
position 30" south (P 3), i.e., 7^ rms — 8.2mK in a 0.16kms * bin 
(Fig.Q. 

3.2. C 18 and 13 C 18 

Example spectra in three isotopes of CO (3-2) are shown 
in Fig. [2] toward two positions in the central region of the 
p Oph A core. Further, Fig. [3] shows the inner, high-resolution, 
map of integrated intensity, §T^dv, of the C 18 0(3-2) line. 
Within a range of R.A. offsets +30" to +50", four distinct in- 
tensity peaks are discernable. In Table[T] these are designated 
P 1 through P 4 and their J2000 coordinates are given. The C 18 
line is very narrow, e.g. merely 1.0 km s" 1 (FWHM) at the in- 
conspicuous (0, 0) position. 

Examination of the entire data set for C 18 (3-2) reveals the 
fact that, within the mapped region, maximum emission occurs 
at LSR- velocities +2.7 to +3.7 km s" 1 . This velocity interval cor- 
responds to t hat of the O? 119 GH z emission viz. v^sr ~ +2.5 to 
+3.5 km s" 1 ( Larsson et al.||2007| and Fig. [6] here. pOphA dis- 
plays a complex velocity field and two distinct velocity com- 
ponents can be identified, giving rise to spectral line blend- 




Velocity 



Fig. 2. (/ = 3 - 2) spectra (r m b vs z^lsr) of, from top to bottom, 
CO (black), C 18 (red) and 13 C 18 (blue) toward two positions 
in the p Oph A core (cf. Fig.j3|. For clarity, two of the spectra are 
offset by +5 K and the 13 C 18 spectra have been multiplied by a 
factor of twenty. 



Table 2. 

T13 C 18 . 



Observed positions of 13 C 18 (3-2) and line opacities, 



Offset 


^LSR 


FWHM 


W 12 c 18 o), 


T 13 C 18 


Note 


(arcsec) 


(kms" 1 ) 


(kms" 1 ) 


W 13 C ls O) (K) 






-60, +80 


2.86 ± 0.05 


0.78 + 0.12 


10.5, 0.25 


0.02 




0, +60 


3.13 + 0.06 


0.92 + 0.14 


8.1,0.22 


0.03 




+30, +80 


3.09 + 0.05 


1.11+0.09 


13.0, 0.38 


0.03 


PI 


+30,0 


3.62 + 0.02 


1.02 + 0.05 


16.5, 0.71 


0.04 


P3 


+60, -100 


2.95 ± 0.06 


1.17 + 0.15 


6.1,0.23 


0.04 


2 lines? 



Note to the Table: Error on line ratios is estimated at 15- 20%. 



ing. These components are essentially confined within the LSR- 
velocity bins [+2, +3] and [+3, +4] (in kms -1 ). Fig. [^presents 
a mosaic of the integrated line intensity in 1.0 km s" 1 wide bins. 
A variety of this kind of channel maps demonstrates quite con- 
vincingly that the location of the O2 emitting gas is most likely 
associated with the central core region of p Oph A. 



4. Discussion 

4. 1 . The dense clumps of p Oph A 

The C 18 intensity maxima in Fig.|3]seem comparable in size 
with the APEX beam, which could indicate that the diameter 
of these clumps does not exceed 20". From the comparison of 
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[+1.0, +2.0] 



[+2.0, +3.0] 




[+3.0, +4.0] 

+ + + :» + v+ + + + 

!+ + + ♦+ Vj^jr + 

:q:: " 



[+4.0, +5.0] 



: 



20 



15 6 

- 10 ^ 




RA offset |"'l 



RA offset f"l 



Fig. 4. A mosaic of maps of integrated C 18 (3-2) line intensity over 1.0 km s" 1 wide velocity intervals, from +1.0 to +5.0 km s" 1 . 
The extended maps with 20" spacing are shown. Two velocity components are identified in the p Oph A core, falling into the [2, 3] 
and [3, 4] bins, respectively. The lowest contour level corresponds to 4 Kkm s" 1 and increments are also by this amount. 



ClB0(3-2) 




RA offset ["] 

Fig. 3. C 18 0(3-2) integrated intensity, fT^dv, of the dark core 
p Oph A. The map was obtained with APEX and observed posi- 
tions are shown as crosses. The beam size at 329 GHz is shown 
in the lower right corner. Offsets are with respect to the origin, 
RA=16 h 26 m 24 s 6 and Dec=-24°23'54" (J2000). The position of 
the outflow driving Class source VLA 1623 is shown by the 
star symbol. P 1-P4 designate the clumps discussed in the text 
(Table[T]). The beam size at 234 GHz is indicated by the dotted 
circles, at the observed 18 positions. 



their loc ations with those o bserved in the emission of the dust at 
1.3 mm ( [Motte et al.fl998) and 850 //m ( [Johnstone et al.|2000| ) 



and of th e quiescent gas in the N2H + (1-0) line ( |Di Francesco et 
al.|2004 ), it becomes evident that P 4 lacks correspondence with 
features at 1.3 mm and N2H + emission, but shows up weakly at 
850yum. P 1 likely is Nl (which is not seen in the dust maps), 
P2 corresponds to N5 (also prominent in the dust as SM IN), 
a nd P 3 seems associate d with N4 and SM 1 (also 16264-2423 
of |Johnstone e t al. 2000 ). Derived temperatures and densities for 
these clumps are of the order of 15-30 K and 0.2 - 5 x 10 6 cm" 3 , 
respectively (e.g., Andre et al.| 19 93 ; Motte et al. 1998 ; Johnstone 
et al.|2000| ). 

In summary, the evidence points toward the fact that also 
O2 is concentrated in the dense dark core regions, where the 
molecules would be protected against photo-dissociation due to 
the intense UV field (Go of the order of 10 2 ) generat ed by the 
two B-stars, east and west of the cores, respectively (Liseau et 
al.|19 99 ). The size of the O2 emitting regions appears not to ex- 
ceed one arcminute, so that a conservative estimate of the Odin 
beam filling would be about 0.01. If the emission originates in a 
core of size ~ 20" or smaller, the Odin beam filling factor would 
be reduced by yet another order of magnitude. The O2 abun- 
dance would scale accordingly and could in this case be locally 
as high as a few times 10" 5 , which would be compara ble to the 
total abundance of oxygen in the gaseous phase (e.g., Liseau & 
|Justtanont|2009] ). 

4.1 .1 . Line optical depths 

The ratio of the 13 C 18 and 12 C 18 line intensities can be used 
to estimate the optical depth in the rarer isotope line, Ti3 C i8 ~ 
ln(l - r l3 )~\ where r 13 = T( 13 C 18 0)/r( 12 C 18 0). From the data 
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presented in Table[2j it is clear that the C 18 0(3-2) line could 
have significant opacity along several lines of sight, unless the 
relative abundance [ 12 C 18 0/ 13 C 18 0] «c 50 (or the excitation 
temperatu res for these spec ies differ substantially). 

Federman et al.| ( |2003| ) determined a column density ratio 
Af( lz CO)/AT( iJ CO) = 125 ± 23 toward a line of sight desig- 
nated pOphA by therrj^] However, their coordinates refer to 
the star, one degree north-north- west from the p Oph A core dis- 
cussed in our paper. In the associated nebula, the physical condi- 
tions are different from those in the dense core, possibly leading 
to different isotopic abundances. In the shielded regions of the 
dense cores, chemical isotopic fractionation can be expected to 
be of minor importance. It is worth noting that in the nearby 
p Oph core C, a lower isotopic ratio [ 12 C]/[ 13 C] = 65 + 10 has 
been derived by |Bensch et al.| ( [2001j ). 



4.2. Column densities 

If local thermal equilibrium (LTE) is a good approximation for 
the level populations, the column density of all molecules of the 
species, Af(mol) in cm -2 , can be estimated from the observed 
intensity of an optically thin line, viz. 



I T l dv 

Af(mol) = — — x 0(r k ) 



fb Imb 



(1) 



where 



o(r k ) = — — — — — —— — —— exp(r u /r k ) (2) 

\ he ) A u i F(T k ) - F(7bg) g u 

with cgs-units of K" 1 cm -3 s. Here, fi> is the beam filling fac- 
tor for the source which may be smaller than the beam (0 < 
fb < 1) and is the main beam efficiency. T tY = hv/k is the 
transition temperature, 7b g = 2.725 K is the temperature of the 
back-ground radiation field, F(T) = T tY /[exp(T tY /T) - 1] is the 
quasi-Planck function, T u is the upper level energy in K, Q(T k ) is 
the partition function and g u = (2/+ 1) is the statistical weight of 
the upper level and the other symbols have their usual meaning. 



4.2.1 . C 18 and H 2 column densities 

We limit the discussion to the central core region, where ob- 
served C 18 0(3-2) line intensities of the +3kms _1 component 
are fT^dv= 20 K km s" 1 . The upper level energy lies nearly 32 K 
above ground. The spontaneous transition probability for the 
transition is A32 = 2.158 x 10" 6 s" 1 , the transition temperature 
is r tr = 15.813 K and the statistical weight of the u pper level 
is g u = 7. Using the collisional rate coefficients of |Schinke et 
al. (1985 ) for collisions with para-H 2 , yields critical densities, 
n C nt ~ A 32 /732Crk), of about 2xl0 5 to 3xl0 4 cm" 3 for T k = 10K 
to 300 K, respectively (Table[3]). Therefore, except perhaps for 
the very lowest temperatures, the condition of LTE should be 
fulfilled for the C 18 (3-2) transition (cf. Sect. 4.1). 

The sizes of the clumps are comparable to the beam size, 
so that the main beam efficiency, rj mh = 0.73, is used for the 
intensity calibration and we assume here a beam filling factor of 
unity. For the broad range of temperatures of 10 to 300 K, the 
corresponding column densities of C 18 are listed in Table^ 
For an X(C 18 0)= 1 to 2 x 10" 7 , the derived H2 column densities, 



on the 20" scale (2400 AU), are N(U 2 ) = (1± 5 05 ) x 10 23 cm 



These results are in general agreement with those reported by 
others ( |Loren et aT1|1990[ |Motte et aL]|1998| ). Possible opacity 
corrections to the C i6 0(3-2) intensity, of the order of r/(l - 
exp -t) > 2, would increase the column density accordingly. The 
column densities presented in Table[3] are therefore likely lower 
limits. 



4.2.2. 18 column density and 18 abundance 

The Nj = 2\ - 0i transition has the largest Einstein coefficient 
of the low-lying 18 Q tr ansitions, viz. A 2 i,oi = 1.33 x 10" 8 s" 1 
(Marechal et al. 1997b). We adopt the coefficients for colli- 
sional de-excitation, 721,01 (7k)> which are based on the work by 



Bergman ( 1995 ) and which have been derived for collisions with 
He. For collisions with H2, these were multiplied by 1.4. Values 
for temperatures other than 300 K were obtained by scaling with 
the square root of the temperature. From Table[3] it can be seen 
that critical densities for the 2\ - 0i transition are rather low for 
a wide range of temperatures (< 1500 cm -3 above 10 K). In par- 
ticular, for the dense core conditions of p Oph A, where densities 
are in excess of 10 5 cm" 3 (Sect. 4.1), LTE is certainly a valid as - 
sumption (see also |Black & Smith|1984l|Marechal et al.|1997b| ). 
The temperature of the transition is T u = T lv = 1 1.228 K and the 
statistical weight of the upper level is g u = 3. 

In Table|3j the results for C 18 and 18 are compared. The 
ratio N(C 18 0)/Af(0 18 0) exceeds unity and increases with de- 
creasing temperature. This ratio could correspond to about half 
the value of that of the CO/0 2 ratio ([Black & Smith|1984| ). For 
three cold cores (10 and 15 K), Fuent e et al.| ( |1993| ) determined 
CO/O2 > 3 - 7, limits consistent with, but considerably smaller, 
than the values displayed in Table[3] The effects of a varying 
C/O ratio in the ISM at column densities (values of the visual 
extinction Ay) as high as those fo und in p Oph A were exp lic- 
itly considered in the models by Mare chal et al.| ([1997a, see 
their Fig. 10). For the O2 1 19 GHz line, the integrated intensity is 
> 10 K km s" 1 for C/O < 0.4 when A v > 20 mag. In contrast, for 
similar extinction, the intensity is < 100 mK km s" 1 for C/O > 1. 
Future observations will likely be able to follow any variation 
of this ratio in diff erent regions of the ISM (see below and also 
|Black&Smith|1984| ). 



Table 3. Column densities of C 18 and 18 





^3,2/73,2 


7V(C 18 0) 


^21,01/721,01 


Af(0 18 0) 


N(C li 


O) 


(K) 


(cm- 3 ) 


(cm- 2 ) 


(cm- 3 )" 


(cm" 2 / 


N(O l 


*0) 


300 


3.4 x 10 4 


4.9 x 10 16 


2.7 x 10 2 


<5.6x 10 15 


>9 




100 


4.0 x 10 4 


2.0 x 10 16 


4.7 x 10 2 


< 2.0 x 10 15 


> 10 




40 


1.7 x 10 5 


1.2 x 10 16 


7.4 x 10 2 


< 9.5 x 10 14 


> 14 




20 


1.7 x 10 5 


1.4 x 10 16 


1.0 x 10 3 


< 6.3 x 10 14 


>22 




10 


1.7 x 10 5 


3.2 x 10 16 


1.5 x 10 3 


< 5.5 x 10 14 


>59 





5 In addition, for their p Oph A line of sight, 
also give Af( 12 C 16 0)/Af( 12 C 18 0) = 1100 ± 600. 



Federman et al. 



(2003 



Notes to the Table: 

a Collision rate coefficient y(T) = 4.9 x 10" 11 (7/300 K)° 5 cm 3 s 
b Limits on Af(0 18 0) are lcr; f h = 1.0 and fT* A dv < 3.3 mK km s" 



In the dense cold (« 100 K) regions of the p Oph A core, 
the column density of 18 is lower than 10 15 cm -2 (Tablepj) 
and, hence, the abundance relative to H 2 , X(0 18 0)< 10 . 
Consequently, for the range of 10 to 40 K and a standard ele- 
mental isotopic ratio, the abundance of the pri mary spe cies of 
molecular oxygen, X(0 2 ) ~ 500/2 x X(0 18 0) ( |Wannier||l980t , 
should be limited to < 5x 10" 7 -2x 10" 6 , consistent with the Odin 
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result (Larss on et al.|2 007). We can conclude, therefore, that in 
the p Oph A core, the molecular oxygen abundance is bounded 
by 5 x 10" 8 < X(0 2 ) < 2.5 x 10" 6 , where the beam averaged 
O2 column density is 10 15 cm -2 . If reflecting the fraction of the 
Odin bea irjj t hat is filled by the O2 source, its implied size is 
<600"/ V250 = 38". This could be well-matched to the 3.5 m 
telescope of the Herschel Space Observatory]^] the beam widths 
of which are 44" at 487 GHz, the frequency of the 2 (Nj = 3 2 - 
h) transition, and 28" at 773 GHz for the (5 4 - 3 4 ) line (Fig.g). 



0.1 



0.01 




0.4 

/(487)//(119) 



Fig. 5. 2 line ratio diagram for the two strongest transitions ac- 
cessible to HIFI aboard the Herschel Space Observatory. In these 
multi-transition calculations, radiation from dust was included 
and LTE was not assumed. Intensity ratios are relative to the (li - 
1 ) 119 GHz li ne, which was detect ed by Odin (24+4 mKkm s" 1 
in a 10' beam; [Larsson et alT|2 007 ). These transitions are the (33 
- 1 2 ) 487 GHz line (44") and the (5 4 - 3 4 ) 773 GHz line (28"), 
respectively. Labels along the graph refer to the density, ft(H 2 ) in 
cm -3 , and the gas temperature, 7\ in K. 



4.3. Nature, location and extent of the 2 source 
4.3.1 . Oxygen in the cold ISM 

The capital l etter d esignation of the cores was introduced by 
Loren et al. ( 1990| ) for the location of emission peaks in lines 
of DCO + in the p Oph cloud. Depending on the details of the 
considered models of the deuteration process, they derived ki- 
netic gas temperatures inside the cores which were always low, 
in the range 18-23 K, whereas temperatures in the outer layers 
were considerably higher. 



6 Odin observations resulted in a column density of oxygen N(0 2 ) = 
1 x 10 15 cm" 2 ( [Larsson et all2007| ). If Q(0 2 ) = Q(C 18 0), the beam 



filling of the 2 source is lO^to 10 2 , i.e. the beam corrected N(0 2 ) ■ 
10 17 - 10 18 cm" 2 . If Af(0 18 0)/Af(0 2 ) = 1/250, then the expected column 
density of isotopic oxygen is likely within 0.4 X 10 15 cm -2 < Af(0 18 0) 
< 4 x 10 15 cm -2 for a source of size about 20" to 60". 
7 |http : //hersche l . esac . esa . int/| 




Fig. 6. The line profile of the C 18 O (3-2) map, after convolution 
with a 600" beam, is shown in black and is compared to the 
scaled O2 line of |Larsson et al. ( 2007| ) shown in red. The higher 
S/N APEX data are spectrally sampled at a higher rate and have 
also higher resolution. The C 18 line intensity is dominated by 
an extended cloud component at an LSR- velocity seemingly dif- 
ferent from that of the P 3 core and the O2 emission (see also 
Tabled 



In the interstellar medium, it is expected that most of the 
molecular oxygen is formed by the reaction 

+ OH^0 2 + H . 

The thermal rate coefficient of this reaction has recently been 
measured in the laboratory down to temperatures of 39 K (Carty 
et al. 2006), where it remains rapid. In subsequent ab initio the- 
oretical calculations, [Xu et aT] ( |2007| ) found a much smaller rate 
at temperatures below 30 K and suggested that this might solve 
the problem of missing O2 in cold interstellar clouds. The low- 
temperature beh aviour of the reaction i s also of interest for ultra- 
cold collisions. Quemen er et al.| (|2009 ) have determined that the 
reaction still proceeds in the limit of zero temperatu re with a 
rate coefficient of approximately 6 x 10" 12 cm 3 s" 1 



Quan et al. 



( 2008) re-examined the sensitivity of the interstellar O2 abun- 
dance to the low-temperature behaviour of the source reaction. 

Also a widely favoured explanation for the generally ob- 
served paucity of molecular oxygen in the gas phase is depletion 
of atomic oxygen with subsequent hydrogenation on cold grain 
surfaces. This scenario seems supported by several observed 
molecules. For instance, the hydrogenation of CO is predicted 
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Table 4. Gaussian parameters of 13 C 18 0(3-2) and 2 (li-lo) 
lines. 



Molecular 


^LSR 


FWHM 




Note 


spectral line 


(kms" 1 ) 


(kms" 1 ) 


(K) 




2 119 GHz 


3.5 


1.5° 


17.4 x lO" 3 


observed 10' beam 


C 18 329 GHz 


3.2 


1.4 


2.9 


convolved 10' beam 


C 18 Q 329 GHz 


3.6 


1.0 


16.5 


observed 19" beam, P3 



a The 2 line is artifically broadened (see Lar sson et al.|2007| ). 

to lead to H 2 CO, CH 3 OH and subsequently also the d euterated 
forms of these species ( Mata r~e7 al. |20 08 ; Fuch s et al.|2009| ). In 
order to become observable, these species have to be returned 
into the gas phase. Indeed, enhanced emission in methanol and 
doubly deuterated formaldehyde has been observed towar d the 
centre of pOphA by, respectively, |Liseau et al. ( |2003| ) and 
Bergman et al. (in preparation). In addition, widespread emis- 
sion of gas phase H2O in pOphA is also observed (Larsson et 
al., in preparation), a fraction of which may have been similarly 
produced by the hy drogenation of O2 on cold grain surfaces 
(Iop polo et al.||20Q8 ). The equilibrium between adsorption and 
desorption of molecules would then naturally lead to low levels 
of both species (as compared and in contrast to pure gas phase 
chemistry). 



4.3.2. Site and size of the 2 source 

Figure[6] shows the C 18 0(3-2) spectrum after the convolu- 
tion of the observed map with a 10' beam. A Gaussian 
profile provides a good fit to the observed line, the pa- 
rameters of which are Tl =2.9 K, ^LSR=3.20kms _1 and 
A^fwhm=1-45 kms -1 (Tableffl). This velocity is offset from that 
of the core P3 (3.62 kms" , cf. also Table[2]) and the inten- 
sity is dominated by an extended component. The integrated 
va lue is fr*di;= 4.5 Kkms" 1 . Using the temperature assumed 
by |Lars son et aL]( |2007| i.e. 30 K) we obtain the beam averaged 
column densit)0Af(C 16 O) > 2.5 x 10 15 cm" 2 . The comparison 
with the Odin result, i.e. N(02) = 1 x 10 15 cm -2 , would indicate 
that N(C n O)/N(0 2 ) > 1. A C 18 abundance that is larger than 
that of O2 would be difficult to explain and would speak against 
an extended O2 emission region. 

The O2 119 GHz Odin line shares the LSR- velocity with 
that of intensity maxima in C 18 and N2H + (this paper and Di 
jFrancesco et al.||2004[ |20Q9| ). It seems therefore reasonable to 
identify the location of predominant O2 emission with the cen- 
tral parts of the cold core p Oph A, i.e. the region including P 2 
and P 3 (SM IN and SM 1, respectively) with a probable extent 
on the 30" to V scale. 



5. Conclusions 

Summarising, we briefly conclude the following: 

• C 18 0(3-2) mapping observations with APEX of a 10' x 5' 
region in p Oph A have revealed a complex radial velocity 



8 This C 18 column density, which represents an average over ten 
arcminutes, implies an H 2 column density, Af(H 2 ) ~ 2. 5 X 10 22 cm" 2 , 
a valu e which has been derived also by other means (La rsson et al.| 
2007 and references therein). At the adopted distance of 120 pc, this 
translates into an H 2 -mass of the p Oph A cloud of > 30 M Q . Not totally 
une xpected, most of the mass would be con tributed on larger scales (cf., 
e.g., |Motte et al.|1998|[Maruta et al.|2009] >. 



field. The central 200" x 200" have been spatially sampled 
at the Nyquist frequency. 

• The z^lsr of the O2 1 19 GHz line appears confined to a par- 
ticular region (SM 1), which is also a prime emitter in C 18 
andN 2 H + . 

• The observation of 18 toward SM 1 (P 3) and SM IN (P 2) 
resulted in upper limits. Combined with the C 18 data, this 
leads to a ratio of Af(C 18 0) to Af(0 18 0) much larger than 
unity. 

• From the O2 and 18 observations we infer an O2 abun- 
dance 5 x 10" 7 < 2 < 2.5 x 10" 6 . 

• The O2 source is likely relatively compact, within a factor of 
a few of 1 arcminute, and should become readily detectable 
by upcoming Herschel HIFI observations. 



Acknowledgements. We wish to thank Cathy Horellou and Daniel Johansson for 
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